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ABSTRACT 

We consider the thermal and non-thermal emission from the inner 200 pc of the Galaxy. 
The radiation from this almost star-burst-like region is ultimately driven dominantly 
by on-going massive star formation. We show that this region's radio continuum (RC) 
emission is in relative deficit with respect to the expectation afforded by the Far- 
infrared-Radio Continuum Correlation (FRC). Likewise we show that the region's 
7-ray emission falls short of that expected given its star formation and resultant 
supernova rates. These facts are compellingly explained by positing that a power- 
ful (400-1200 km/s) wind is launched from the region. This wind probably plays a 
number of important roles including advecting positrons into the Galactic bulge thus 
explaining the observed ~kpc extension of the 511 keV positron annihilation signal 
around the GC. We also show that the large-scale GC magnetic field falls in the range 
^100-300 /iG and that - in the time they remain in the region - GC cosmic rays do 
not penetrate into the region's densest molecular material. 

Key words: cosmic rays - ISM: jets and outflows - supernova remnants - Galaxy: 
centre - radio continuum: ISM - gamma-rays: theory 



1 INTRODUCTION 

The extreme ISM conditions in the central ~ 200 pc of the 
Galaxy render the region more akin to a star-bursting sys- 
tem (e.g., |Launhardt et aL] [2002) than to almost any re- 
gion in the Galactic disk. The similarities include: i) a high 
areal star-formation and (consequent) supernova rates; ii) a 
flatish overall radio spectrum within the star-forming region 
(cf. |Niklas et aL]|1997| |Thompson et al.||2006[ ); iii) a region 
surrounding the star-forming nucleus of bright but diffuse, 
non-thermal radio emission; iv) the existence of diffuse 7-ray 
emission also apparently associated with star formation (cf. 
Fermi LAT Collaboration||2009l |HESS Collaboration|[20u9l 
VERITAS collaboration||2009| on NGC 253 and M82); and 
v) a rather strong magnetic field (> 50 pG; Crocker et al. 



20101. Here we argue for another similarity: a strong out- 



flow with a speed 400-1200 km/s (comparable to the es- 
cape speed) and energetically consistent with being driven 
by current star- formation ( Veilleux et al.|200 5 Strickland & 
|Heckman||2009|). 
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Massive, young stars are copious producers of UV 
and optical light which is reprocessed into IR emission 
by the dust of the stars' natal molecular envelopes (|De-| 



vereux & Young 19901. On the other hand, cosmic ray 
(CR) electrons and ions - ultimately powered by super- 
novae (e.g., Hillas 2005) - produce their own (non-thermal) 
radiative signatures. These include ~GHz radio contin- 
uum (RC) synchrotron emission and inverse Compton (IC) 
and bremsstrahlung emission at 7-ray wavelengths by CR 
electrons and 7-rays from neutral meson decay following 
hadronic collisions between CR ions and gas. 

Given the connection of these radiative processes back 
to massive (A/* > 8M Q ) star formation ( |Volk||1989| ), one 
might expect that they be globally correlated. Such is ob- 
served (Dickey & Salpeter 1984; d e Jong et al.|1985 Helou 



et al.|[i"985| ): an extremely tight (dispersion of ~ 0.26 dex 
Yun et al. 200 1) FIR-RC correlation (FRC) is found (e.g., 
Condon||i"992"t to hold over five orders of magnitude in RC 



luminosity ( |Yun et al.|2001[), and b oth globally and at sub- 



galactic scaies (H ughes et al.|2006 



Tabatabaei et al.|2007l 



2006 



Likewise, one might also expect (Thompson et al 
|Thompson et aLj|2007| a global scaling between FIR and 
7-ray production ('F7S'). As we show below, however - and 
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in interesting contrast to star-bursting systems (Thompson 
et al.|2006 1 - the GC does not fall on these scaling relations: 
we detect far less non-thermal emission than expected given 
the region's star-formation rate. This deficit is ultimately 
explained by a large-scale, powerful outflow from the region. 



2 CORRELATIONS AND SCALINGS 

The H.E.S.S. Imaging Air-Cherenkov 7-ray Telescope has 
detected hard-spectrum, diffuse ~TeV 7-ray emission sur- 
rounding the GC over the region defined by |Z| < 0.8° and 
|6| < 0.3° with an intensity of 1.4 x 10~ 20 cm -2 eV" 1 s _1 
sr _1 at 1 TeV (with the point TeV source coincident with 
Sgr A* subtracted). Only dimmer diffuse TeV emission is 
detected outside this (hereinafter) 'HESS field'. 

Of note, is that the spectral index, 7, of the GC diffuse 
~TeV emission, where F 1 oc E' 1 , is 2.3 ± 0.07 5tat ± 0.20 sys , 
significantly harder than the spectral index of the CR ion 
population threading the Galactic disk and the diffuse 7- 
ray emission it generates. Disk CRs experience energy- 
dependent confinement and their steady-state distribution 
is, therefore, steepened from the injection distribution into 
the softer ~ E~ 2 ' 75 spectrum observed at earth (see, e.g., 
I Aharonian et~aT||2006) > . The GC TeV 7-ray spectral index 
(and that inferred for the parent CR ions) is close to that 
inferred for the injection spectrum of Galactic disk CRs, 
itself within the reasonable range of ~ 2.1-2.2 expected 
(Hillas 2005) for lst-order Fermi acceleration at astrophysi- 
cal shocks. 

Empirically the 1.4 GHz RC (spectral) luminosity and 
the total IR luminosity (Ltir[8 — 1000] \im\ Calzetti et al. 
2000 1 are connected as ( Yun et ah||2001| |Thompson et alT 
20071 



vL v (lA GHz) ~ 1.1 x 10~ 6 Ltir 



(1) 



with a scatter of ~ 0.26 dex. On the basis of IRAS data 
(Launhardt et al.l |2002h the Ltir of the HESS field is 



1.6 x 10*" erg/s, implying ( Kennicutt 1998) a SFR of 0.08 



M /yr for the HESS field. In useful units, the 1.4 GHz RC 
luminosity ( |Reich et al.|l990| l of the HESS field is 1.7 x 10 35 
erg/sQ, dex or ~ 4<r short of the expectation from the 

FRC. 

|Thompson et al.|p007 l use the empirically-established 
connection between the SFR and the total infrared luminos- 
ity to relate the power, injected by supernovae into CRs, to 
Ltir and hence to predict that the TIR and 7-ray emission 
from luminous star-forming galaxies should scale as 



^L„(GeV) ~ 2.0 x 10 5 7/0.10 Ltir 

where the proton spectrum is assumed oc Ep 2 up to E. 



(2) 



10 lb eV and we have renormalized the equation of [Thomp- 



son et al. (20071 assuming 770.10 10% of the 10 51 ergs per 



supernova goes into relativistic ions. This relation assumes 
that the region under consideration is calorimetric to CR 
ions. 



On the basis of the results presented by Meurer ( 2009 1 , 



1 We have removed the contribution from synchrotron emission 
from relativistic electrons in the Galactic plane but out of the 
GC: see Crocker, Jones, Aharonian et al. 2010, to be submitted, 
henceforth Paper II. 



Fermi observes a luminosity of ~ 3 x 10 36 erg/s for E 1 > 
GeV for emission from the central 1° x 1° field, only ~10% 
of that expectated from the FIR emission. The Fermi obser- 
vations are, however, substantially polluted by line-of-sight 
and point source emission (including from a source coinci- 
dent with Sgr A*: Chernyakova et al. 20101 so they only 



constitute an upper limit to the true diffuse 7-ray emission 
from the region. 

We can consider the HESS data by scaling eq. [2] from 
L 1 {E 1 > GeV) to L~,{E 1 > TeV). For the TeV spectral in- 
dex of ~ 2.3 and assuming an hadronic origin to the TeV 
7-rays, L 1 {E 1 > TeV) ~ 0.2 L 1 (E 1 > GeV). The TeV lu- 
minosity we infer for the HESS field of 1.2x 10 35 erg/s (inte- 
grating to 100 TeV) is only ~£% of the prediction from the 
suitably-scaled version of eq. [1[ 

Thus the FRC fails badly in the case of the HESS field: 
far less RC than expected is detected given its FIR output. 
Likewise, the 7-ray luminosity of the region is significantly 
in deficit given the region's FIR output (and implied star- 
formation rate). There are three potential explanations of 
these discrepancies: 

Firstly, a RC deficit could arise if a starburst event oc- 
curred more recently ( £ 10 7 years) than the lifetime of the 
massive stars whose supernova remnants accelerate the CR 
electrons which generate synchrotron emission. Although we 
expect some stochastic variation in the GC's overall SFR, 
we find, however, that the current SFR is close to the long- 



term ( > 10 7 ) average value (cf. |Serabyn fc Morris 



1996 



Figer et al.|2004 1 . A strong piece of evidence for this is that 



a number of other handles on the GC supernova rate we 
describe in Paper II that are sensitive to long-term aver- 
age values of this quantity - through, e.g., studies of the 
region's pulsar population ( Lazio fc Cordes|2008 1 - are con- 
sistent with the supernova rate implied by the current SFR 
as traced by FIR, viz. 0.04/century in the HESS field. 

Secondly, it may be that GC SNRs are intrinsically low- 



efficiency (cf. Erlykin & Wolfendale 20071 CR accelerators 
(plausible because of their - on average - dense environs: 
Fatuzzo & Melia 20051. However, the detailed numerical 



modelling set out in Paper II shows that GC supernovae 
do, indeed, accelerate CRs with typical (e.g., |Hillas|[2~005[ ) 
efficiency: about 10% of the total 10 erg mechanical en- 
ergy per supernova goes into non-thermal particles. Given 
the above rate, this implies that supernovae inject ~ 10 39 



erg/s into the GC CR population (cf. Crocker & Aharonian 



2010) 



Lastly, given the half-height of the region is only ~ 40 
pc, a reasonable reaction to the breakdown of the FRC is 
that it is unsurprising; many studies ( Murgia et al. 2005 1 



find a break-down in the correlation at ~ kpc, often pro- 
posed to be due to electron transport. On the other hand, 
studies, e.g., of the Large Magellanic Cloud (Hughes et al. 



20061, the Scd galaxy M33 (Tabatabaei et al. 20071, and 



within the Milky Way (Zhang et al. 20101 reveal a tight 



connection between RC and FIR emission down to scales 
< 50 pc. 

A potential fourth explanation of why the HESS field 
falls off the FRC is that power fed into non-thermal elec- 
trons is 'lost' to ionization and bremsstrahlung and/or in- 
verse Compton emission ( Thompson et al.|2006 Thompson 
|et al.|2 007) rather than synchrotron emission (plausible be- 
cause of the GC's dense gas and radiation environment) 
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Given, however, the HESS field also falls short of the F7S, 
this explanation is, at least, seriously incomplete. 

In summary here, it seems that CR transport out of the 
HESS field is by far the most plausible explanation for why 
it falls off the global scalings discussed; below we show that 
the transport mechanism is a wind. 



3 PRIOR EVIDENCE FOR AN OUTFLOW 
FROM THE GC 

There is multi-wavelength evidence in support of the ex- 
istence of GC outflow. Recent infra-red observations show 
that the GC's massive stellar clusters are blowing a bubble 



into their environment (e.g., Bally et al. 20101. Keeney et 



al. ( 2006 1 and Zech et al. ( 2008 1 have found evidence for 



high-velocity gas consistent with a GC outflow or fountain 
in UV absorption features towards, respectively, two AGN 
and a GC globular cluster. The region's spectacular non- 



thermal radio filaments (Yusef-Zadeh et al 
due to a fast outflow (e.g., Shore & LaRosa 



1987] ) may be 
1999T. RC evi- 



dence of an outflow was found in 10 GHz radio continuum 
emission by |Sofue fc Handa ( 1984 1 in the form of a ~ 1° 
(or ~ 140 pc) tall and diameter < 130 pc shell of emission 
rising north of the Galactic plane called the Galactic Centre 
lobe (GCL). RC emission from the lobe's eastern part has 
HI absorption that clearly puts it in the GC region ( Lasenby 
et al.|1989 1 and its ionized gas has a high metalicity ( Law et 
al.||2009^ ~Filamentary structures coincident with the radio 



have been discovered at mid-infrared wavelengths (Bland 



Hawthorn & Cohen||2003 1 and the structure interpreted as 



evidence for a previous episode of either st arburst (|Bland- 
Hawthorn & Cohen 2003 ) or nuclear activity ( Melia & Falcke 



2001). Law (20101 has found that the formation of the GCL 



is consistent with currently-observed pressures and rates of 
star-formation in the central few x 10 pc of the Galaxy. Fi- 
nally, Law (20101 determined the ~GHz spectral index of 



the GCL steepens with increasing distance (both north and 
south) of the Galactic plane. This constitutes strong evi- 
dence for synchrotron ageing of a CR electron population 
transported out of the plane. Thus, a natural interpretation 
is that the GCL's RC emission is due to CR electrons ad- 
vected from the inner GC (essentially the HESS region) on 
a wind (cf. |Zirakashvili fc Volk|2006[ |Heesen et al.|2009] on, 
e.g., NGC 253 and M82). 

This interpretation requires that 

• The spectrum of the electrons leaving the HESS region 
(as given by Eq[3| must match the spectrum at injection 
required for the GCL electrons with spectral index 2.0 - 
2.4 (ICrocker et al.||2010h. This will be well-satisfied if an 



energy-independent transport process like a wind removes 
CR electrons - accelerated into an in-situ ~ E~ 2 distribu- 
tion - from the inner GC. 

• The power in electrons leaving the HESS region must 
be enough to sup port the GCL electro n population, viz. (3 — 
10) x 10 37 erg/s ( |Crocker et aL|2010| ). This is well satisfied 
given the SN rate in the HESS region. 

• The time to transport electrons over the extent of the 
GCL must be less than the loss time over the same scale. 
This implies a wind speed of strictly > 150 km/s and prob- 
ably > 300 km/s: see fig. [t] 
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Figure 1. Lower bound on the wind speed required for electrons 
advected out of HESS field to synchrotron-illuminate the entire 
extent of the GC lobe within their loss times given by ionization, 
bremsstrahlung, synchrotron, and IC emission for environmental 
parameters of B and njj and an interstellar radi ation field en- 
ergy d ensity Ujsrf — 20 eV cm -3 . We infer from |Ferriere et al.| 
that the volumetric average njj in the GCL is ~ 10 cm 3 . 
The strict lower limit to the GCL magnetic field at 50 fiG (and 
probable value 100 fiG: |Crocker et al.|2010[ | imply a conservative 
lower limit to the GC outflow speed of > 150 km/s (and probably 
£ 300 km/s). 



4 NON-THERMAL HINTS OF AN OUTFLOW 
FROM THE GC 

An important consideration is why the GC CR ion popula- 
tion is so hard in comparison to the diffusion-steepened, lo- 
cal population. There are three reasonable interpretations of 
this: i) the system is out of steady state with less time having 
passed since the CR injection event than required for diffu- 



sion steepen ing (cf. the interpretation adopted by Aharonian 
et al.|2006 that a single CR-injection event ~ 10 year ago 
at the GC explains the observed diminution in the 7-ray to 
molecular column ratio beyond \l\ ~ 1°); or there is a small- 
est relevant timescale defined by an energy-independent ii) 
loss or iii) escape process. 

We argue here that iii) is preferred by all the evidence. 
We can dismiss ii) on the basis of our results above which 
show the system falls far short of being a calorimeter for 
protons. A number of factors also tell against i): firstly, as 
argued, other evidence indicates that the system is close to 
its steady state; secondly, the spectral index of the ~TeV 
emission is a constant ~ 2.3 over the HESS region (within 
errors) presenting, therefore, no evidence of diffusion hard- 
ening at the leading edge of a (putative) diffusion sphere; 
and lastly, these spectral considerations apply also to the rel- 
ativistic electron population: the hard radio spectrum of the 
region, a 0.54 (for S v oc v~ a and radio data 1.4-10 GHz: 
see Paper II), requires that the synchrotron-emitting elec- 
tron population is also very hard, ~ Eg 2 ' 1 . Given the rather 
short loss times associated with synchrotron and IC emission 
in the GC environment, this hard electron spectrum con- 
stitutes independent evidence for rather quick and energy- 
independent CR transport (cf. Lisenfe ld fc Volk|2000 |. 



Consider then the region's non-thermal particle popu- 
lation which, in steady state, approximates to: 



n x (E x 



Ti oas (E x )T asc 



T loa £E x ) + (7 - l)r oac 



x(E x ) 



(3) 
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where Q e (E e ) denotes the injection rate of particles of type 
x g {e,p}; we account for both escape and energy loss over 
t csc and Ti osa with the escape time assumed to be energy- 
independent; and 7 is the spectral index of the (assumed) 
power-law (in momentum) proton or electron spectrum at 
injection. 

Turning now to the CR ion population (henceforth pro- 
tons for simplicity), we have already seen that we only detect 
~2% of the TeV 7-ray flux expected in the calorimetric limit. 
Given, then, that pp collisions are by far the dominant en- 
ergy loss process for high-energy CR protons, this deficit im- 
plies that there is significant escape of accelerated ions (with 
accompanying adiabatic losses) - i.e., the system is quite 
far from calorimetric. We define R T cv = L^v/L^w* ~ 10~ 2 
(uncertain by a factor ~2) as the ratio of the observed flux 
of TeV 7-ray emission to the expected in the calorimetric 
limit (cf. fractions ~ 0.01 and ~ 0.05 for the Galactic disk 
and NGC 253; |HESS Collaboration||2009} . From Eqfjand 
accounting for adiabatic losses with timescale 7~f dbtc = 3t^ sc : 



10" 



3t„ p 



3t£c + 4r„. 



(4) 



By analogy with the hadronic case, we define R la 



J LI 



10" 



(again uncertain by a factor ~2). 



Given the very flat radio spectral index, this deficit is po- 
tentially explained as a result of electron energy loss into 
bremsstrahlung, adiabatic deceleration or advective escape. 
Eqfjthen gives 

- c e ac (r brcms +3r,f 3C ) 
y »ch(r brcms +4r c e ac ) 



R 



0.1 



(5) 



Now, given the foregoing, particle escape is both energy- 
independent and the same for CR electrons and protons 
(t^ c = t* c = const) as would be expected for a wind. This 
means that eq. [4] and [5] yield a combined constraint on the 
required velocity of the outflow responsible for particle re- 
moval: see fig. [2] 

Also shown in fig. [2] are minimum and maximum val- 
ues for the speed of the star-formation-driven 'super-wind' 
expected on the basis of observations of the nuclei of ex- 
ternal, star-forming galaxies and the GC's high areal SFR 
(Strickland & Heckman 2009). The asymptotic speed of such 
a wind scales as v w i„ d ~ a/2 77 E/M where < 77 < 1 is 
the thermalization efficiency, typically ranging between 0.1 
for relatively quiescent star formation and almost 1 for star- 
bursts (Strickland & Heckman 2009). Adopting r\ m i n = 0.1, 
r] max = 1.0, E = 1.4 x 10 4U erg/s and 0.025 M /year (see 
Paper II), we find v™^ d 400 km/s and w^Td - 1200 km / s - 

Putting some of these considerations in a different form, 
we expect a TeV luminosity from the HESS region which 
satisfies L 1 (E 1 > TeV) ~ 1/3 U C n(E p > 10 TeV)/r pp V < 
L; b5 (£ 7 > TeV) = 1.2 x 10 35 erg/s where U CR {E P > 
10 TeV) ~ 1/20 x 1.4 x 10 39 erg/s x d/v wind /V is the en- 
ergy density in CR protons sufficiently energetic to generate 
TeV 7-rays, d ~ 40 pc and V ~ 10 62 cm -3 for the HESS re- 
gion, and riH is the effective gas density the protons sample. 
This implies nn 6 cm" 3 (v W i n d/1200 km/s), cf. the volu- 
metric average gas density through the HESS region ~120 
cm" 3 summing over all phases and ~6 cm -3 including only 
plasma phases. Likewise, the total gas mass the protons sam- 
ple satisfies M gas 5 x 10 5 Mq (v wind /12QQ km/s) which is 
much less than the ~ 10 7 Mq of gas in the region. In order 
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Figure 2. Outflow speed inferred given the departures from 
calorimetry for both protons ('p') and electrons ('e'): -Rtcv = 0.01 
and -Rradio — 0.1 as described in the text (the width of the bands 
reflects the uncertainty of ~2 in both i?Tev and i? ra dio). Protons 
cool via their hadronic collisions with ambient gas (hence the lin- 
ear dependence between wind speed and gas density, nu) and 
adiabatic deceleration. In addition to bremssrahlung (and ioniza- 
tion) , electrons also cool via synchrotron (so the magnetic field en- 
ters as a parameter) and IC emission and adiabatic deceleration. 
As the wind escape time is the same for both electrons and pro- 
tons, the intersection of the electron and proton bands describes 
a valid gas density and wind velocity for the HESS environment 
for each magnetic field sampled. The horizontal dashed line shows 
the approximate maximum allowed wind speed (~ 1200 km/s) 
balancing the total power assumed injected into the system by 
supernovae and massive stars (1.4 X 10 40 erg/s) with the kinetic 
power advected by the wind plasma at its asymptotic velocity 
(assuming 100% thermalization efficiency). The horizontal dot- 
dashed line shows the approximate minimum plausible wind speed 
(~ 400 km/s) for thermalization efficiency of 10%. 



that the region's protons not sample all the molecular gas 
in the region they should be removed in a time shorter than 
the convection time into the dense regions of the molecular 
clouds: twind = d/v wind < t c i oud ~ 10 pc/30 km/s (adopting 
30 km/s as a typical internal velocity dispersion for the re- 
gion's giant molecular clouds, conservatively, of radius ~10 
pc: e.g. |Morris fc Se rabyn 19961 which also implies a lower 
limit: v w i nd ^ 130 km/s. Typical timescales are plotted in 

fig-Hi 



5 DISCUSSION AND CONCLUSIONS 

A clear picture emerges from the considerations above. 
Given the morphological and spectral data on the GC lobe, 
we can infer that it is illuminated with CR electrons injected 
in the HESS region carried from the plane on an outflow with 
a speed 150-1000 km/s. The spectral data on the HESS re- 
gion itself imply that most CR electrons and protons accel- 
erated in situ are advected from the region; electrons lose 
only O[10 %] of their power to synchrotron emission in the 
HESS region while protons lose only 0[1 %] of their power 
to pp collisions on ambient gas in the same region. 

Self-consistently and given our understanding of out- 
flows from external, star-forming galaxies, the same star- 
formation and subsequent supernova processes that drive 
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Figure 3. HESS region timescales for central parameter values 
suggested by our analysis, viz. njj = 10 cm -3 , and u w i„d = 700 
km/s with i) (horizontal solid band) the inverse of the super- 
nova rate; ii) (dashed horizontal line) particle escape with energy- 
independent velocity of 700 km/s; iii) (solid red lines) electron 
cooling for (thick) B = 2 X 10 -4 G and (thin) the limiting case 
of vanishing magnetic field (IC cooling dominant at high energy); 
iv) (blue dotted line) proton cooling. Calorimetry generically re- 
quires t loss < t esc . 



the thermal and non-thermal radiation from the HESS re- 
gion will also drive an outflow with a speed 400-1200 km/s. 
This implies that the magnetic field in the HESS field lies 
in the range 100-300 /iG and the effective gas density en- 
countered by the CRs is in the range 3-20 cm -3 . The latter 
is much less than the volumetric average Uh over the HESS 
region suggesting that even super- TeV CRs do not 'sample' 
all H2 before escaping the region. 

We suspect that the outflow we identify plays many 



important roles (see Paper II and Crocker & Aharonian 
2010) including advecting positrons into the Galactic bulge 
(thereby explaining the ~kpc extension of the 511 keV an- 
nihilation radiation: Weidcns pointner et al.|[2008 l, carrying 
CR ions accelerated by GC supernovae out to very large 
heights (~ 10 kpc) thereby explaining the WMAP 'haze' and 
Fermi 'bubbles ' (|Finkbeiner et al-pOOl) [Dobler et al.||2010 



Su et al. 2010 Crocker & Aharonian 20101, and generally 



keeping the energy density of the non-thermal components 
of the GC ISM in check ( |Breitschwerdt et al.|2002 i. 
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